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Classical relativistic system of point particles coupled with an electromag- 
netic field is considered in the three-dimensional representation. The gauge 
freedom connected with the chronometrical invariance of the four-dimensional 
description is reduced by use of the geometrical concept of the forms of rel- 
ativistic dynamics. The remainder gauge degrees of freedom of the electro- 
magnetic potential are analysed within the framework of Dirac's constrained 
Hamiltonian mechanics in the front form of dynamics. The results are im- 
plemented to the problems of relativistic statistical mechanics. Based on the 
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tem is written down in a gauge-invariant manner and an integration over field 
variables is performed. 
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1 Introduction 



All fundamental interactions in physics have got gauge nature and demand the 
use of singular Lagrangians. An adequate method for dealing with such systems 
was developed by Dirac (see [0, |], ^). His approach, i.e., constrained Hamiltonian 
mechanics, has been elaborated in many directions |[ and usefully applied to 
the problems of quantum field theory (see, e.g., 0, But an application of this 
technique to the relativistic statistical mechanics ^ seems rather obscure, al- 
though, for example, the thermo field formulation of the condensed matter theory 
|T0| admits a variational approach with the same gauge structure as in quantum 



electrodynamics and in related theories. We only note the paper |Tl|, where the 
use of Faddeev' measure was suggested for the construction of the relativistic 
partition function of directly interacting particle system described by means of the 
constrained Hamiltonian formalism. 

The present paper is concerned with the classical relativistic system of point 
charges coupled with an electromagnetic field. The first attempt to analyse the 



constraint contents of such a system was made by Dirac The corresponding 



action functional in the four- dimensional representation has got two kinds of the 
gauge freedom. The first is connected with arbitrariness in the parametrisation 
of particle world lines (chronometrical invariance); the second is generated by the 
proper gauge transformations of electromagnetic potentials. We reduce the gauge 
freedom of the first kind by means of the farther great Dirac's invention, namely, the 
concept of the forms of relativistic dynamics [|1^ . Here we use the notion of the form 
of dynamics as denoting the description of relativistic system, which corresponds 
to a given global simultaneity relation defined by means of some foliation of the 
Minkowski space by space-time or isotropic hypersurfaces. Moreover, we consider 
only the case when the corresponding simultaneity relation is independent of particle 
or field configurations (cf. Ijl^)- Using Dirac's constraint formalism, we try to isolate 
the gauge degrees of freedom and to formulate the statistical description of the 
system in a gauge-invariant manner. 

The general features of the formalism were presented in |]15| within the instant 
form of dynamics. It will be our purpose in the present paper to explore the further 
possibilities connected with other forms of relativistic dynamics; in particular, we 
shall consider the front form description of the charged particle system. 

The paper is organised as follows. In Sec. 2 we will first establish the structure 
of Hamiltonian description of the system of particles plus field in an arbitrary form 
of relativistic dynamics. Our results are then applied to the front form of dynam- 
ics. In Sec. 3 we perform the analysis of the corresponding constraints. Section 4 
contains the elimination of the gauge degrees of freedom by suitable canonical trans- 
formation. Sections 5 and 6 are devoted to the application in relativistic statistical 
mechanics. We formulate a Liouville equation for distribution function in the front 
form of dynamics and write down its equilibrium solution corresponding to classical 
Gibbs ensemble accounting particle and field degrees of freedom. Peculiarity of the 
front form of dynamics allow us to perform an integration over field variables in a 
relativistic partition function. Some calculations of a purely technical nature are 
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collected in appendices. 



2 Charged particle system in an arbitrary form 
of dynamics 

We shall consider a system of charged particles, which is described by their (time- 
like) world lines in the Minkowski space-timeQ addtoresetequationsection 

7, : R ^ M4, r ^ x^(r). (2.1) 

An interaction between charges is assumed to be mediated by an electromagnetic 
field F = dA with the electromagnetic potential 1-form 

A = A^{x)dx''; (2.2) 

F = ^Ff,,{x)dx>' A dx\ F^,{x) = d^A, - d,A^, (2.3) 
= d/dx^ . The dynamical properties of such a system are completely determined 



by the following action functional |jT6|, |T^, [18 

S = J2l {-niaV-ulir,) + eaU:{ra)A,[xa{ra)]] - J d^xF,,{x)F^''{x), 

a=l 

(2.4) 

where ma and denote the mass and the charge of particle a, respectively, u^{Ta) = 
dx'^ija) / dTa- We are interested in constructing the Hamiltonian description of the 
system in a given form of relativistic dynamics. The particle and field degrees of 
freedom will be treated on equal level. 

The form of relativistic dynamics in its geometrical definition is specified by the 
foliation E = {Ej | t G R} of the Minkowski space-time with space-like or isotropic 
hypersurfaces 

St = {x G M4 I a{x) = t}; (2.5) 
V'"'id^cr)id,a)<0. (2.6) 



see 



19, 20, 21]). As it follows from condition (|2.6|) , doa 7^ and therefore the 
hypersurface equation (|2.5|) can be solved with respect to in the form: 



a;° = (^(t,x), ^={x\x^,x^). (2.7) 

For definiteness we put 

doCT>0 and ^^^i^ = (^,(t,x) >0. (2.8) 



^The Minkowski space-time M4 is endowed with a metric = diag(— 1, 1, 1, 1). The Greek 

indices /z, i', . . . run from to 3; the Latin indices from the middle of alphabet, fc, . . ., run from 
1 to 3 and both types of indices are subject of the summation convention. The Latin indices from 
the beginning of alphabet, a, b, label the particles and run from 1 to N . The sum over such indices 
is pointed explicitly. The velocity of light and the Planck constant h are equal to unity. 
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Making use of the identities 



•^oO-Ie, = V^t ^ dia\j.^ = ipi = dip/dx\ (2.9) 
it is easy to see that condition ( ^.61 ) imphes the following inequality for the function 

f dip dip dip 



Then one can consider the particle world lines 7^ as being determined by a set 
of points Xa{t) = 7an^t of intersections with the elements of the foliation E. The 
parametrical equation of the world line in a given form of dynamics is 

X° = Xl{t) = p{t,^a{t)) = Pa, = (2.11) 

The variable t serves as a common evolution parameter of the particle system. 
Therefore, the choice of the form of dynamics is equivalent to the fixing of the 
parameters Tq = t of the particle world lines in the reparametrization-invariant 
action ( ^ . 

On the other hand, the foliation S specifies a certain 3 + 1 splitting of the 
Minkowski space-time M4 = M x Sq, because by the definition of a foliation all 
hypersurfaces for different t are diffeomorphic to Sq. That splitting, / : M4 — >■ 
M X So, is really determined by (|2.7|): 



/:(x°,x)h^(^(t,x),x). (2.12) 
Accounting that the determinant of the transformation ( p. 12]) is 

9(x°, x) 



a(t,x) 



pt{t,^), (2.13) 



one can immediately rewrite the action functional (|2.4]) into the form 



S= dtL (2.14) 



with the Lagrangian function 

N 



L = J2 {-maViDpa)^ - v2 + [Ao{t,^a)Dipa + vlA{t,^a)]}+ [ d'xO, (2.15) 

a=l 

'C = --^Mt,^)F,.F^^''. (2.16) 

Here = dxa/dt, D = d/dt = d/dt + Y,^ v\dldx\, and A^ = A^o /"i, F^^ = F^^ o 
f~^. The dynamical variables of our variational problem will be the functions Xa(t), 
yl^(t, x) and their first order derivatives with respect to the evolution parameter, 
Va(t) and y4^^j(t,x). Using the obvious relations 

9oi^ = A^,iV?r\ diA^ = A^,^i - A^^t'^ip-^ , (2.17) 
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we find the following expression for the field Lagrangian density: 

£ = ^ {e,e' - h,h') , (2.18) 

where 

ei = -Ao,i + iA^t + Ao^tVi)ip;\ h,=eijk{A^'^ - A\ip^ip;^) . (2.19) 
The canonical momenta of our problem are given by 

Vai = Tr- = ^ + Ca [Aiit, Xj + Ao(t, ^a)^ai\ , (2.20) 

dC 1 

i7,(t, x) = = - (e. - e,,,^^h^) , (2.21) 

The basic Poisson brackets are 

{A{t, x), 77, (t, y)} = 6,,6'{^ - y), (2.23) 
{Ao(t,x),i7(t,y)} = -53(x-y); 

all other brackets vanish. Equations (|2.21|) and (|2.22|) imply the primary constraint 

^ = n + ilVi = 0. (2.24) 

In the isotropic forms of dynamics, however, some additional primary constraints 
may occur. To see that, we rewrite equation (|2.21| ) into the form 

Hiit, x) = ^ i&i + PijA^t^;^ + Ao^t^i^;^) , (2.25) 

where &i does not contain the derivatives with respect to the evolution parameter 

&i = -Ao^i - ip^hij] hij = Aj^i - Aij, (2.26) 

and 

A, = (1 - ^^ij + Vi^r (2.27) 
Upon using the identity I3ijip^ = ipi, one obtains equation 

A, {A\t + Ao,t^') = ipt (47ri7, - 6^) . (2.28) 

The matrix j3 = has determinant det /? = (1 — (p'^y. 

In the case of space-like forms of dynamics (c^^ < 1) the matrix is nonsingular 
and possesses an inverse 
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so that equation (|2.28|) can be solved as 

A,i + Ao,m = -iij^t (4vri7^^ - &) . (2.30) 

For isotropic forms of dynamics we have (3ij = (piipj and (|2.28| ) imphes constraints 

6 = e^jk^' {Ann'' - e') = 0. (2.31) 

Because it holds identically ip'^C.i = 0, really only two additional primary constraints 
occur. 

Next, we consider the canonical Hamiltonian of our system, which is defined as 



N 

He = Y,P-^< + / {^'^^,t - nAo,t) - L. 

a=l ^ 



(2.32) 



The immediate calculations give 



N 

Hc = J2Ha{t,Xa,Pa)+ / d^xHc 
a=l ^ 



with 



and 



- ea9?at^o(i,Xa) 



(2.33) 



(2.34) 



(2.35) 



The problem of solution of ( 2.20 ) with respect to the particle velocities is analogous 
to the case of free particles [Q. For completeness we collect the corresponding 
expressions in Appendix A. Replacing into ( |A.9|) Pi by 

kai = Pai - ea [Ai{t, X^) + Ao{t, :>Ca)'-Pai] (2.36) 

yields the one-particle Hamiltonian for the space-like forms of dynamics: 



V/(1 - cpDiml + kl) + [kiipaif + KVa] - eaVatAo(t,^a). (2.37) 



Similarly, for the isotropic forms of dynamics we obtain from ( |A.11D 

1 ml + 'kl . . . 

Ha = --^at—Ti ea^atAo{t,Xa). 

2 Kfai 



(2.38) 



To express the field Hamiltonian ( p.35| ) in the terms of canonical variables we 
rewrite it into the form 



He 



1 

8^ 



1 



(2.39) 



6 



Then for the space-hke forms of dynamics, taking into account ( 2.29 ) and (|2.3CI|) , 
we obtain 



= ^Vt (47ri7, - Gi) {Ann, - &,) + h,,h'^ + Ao,iA 



OA 



(2.40) 



In the case of isotropic forms of dynamics (c^^ = 1), when Pij = ipifj, ( |2.28| ) gives 

ip, {A^,t + Ao,tv') = ipt^, {AnW - &) , (2.41) 
so that ( |2.39| ) takes the form 

K = ^V't I [p' (4vri7, - ©i)] ' + hii.h'^ + Ao.AO'* I . (2.42) 

Detailed analysis of the Hamiltonian description for the charged particle system 



in the instant form of dynamics was carried out in [|l^ with the application to 
the classical relativistic statistical mechanics. In the following sections we shall be 
concerned with the front form of dynamics, analyzing additional primary constraints 
and ensuring the corresponding Liouville equation. 



3 Constraint analysis in the front form of dynam- 
ics 

Let us consider the family of the forms of dynamics, which is given by 

x° = t + nx, = 1. 
According to ( ^ ), dOlp , the set of primary constraints is 



(3.1) 



i = n + Wrii = 

= Eijkn' {Ann' + Ao'' - riih'') = 0. 



It is easy to see that ( |3.3| ) is equivalent to 

{Ann'' + Ao''')_ 



Ujh^'' = 0, 



(3.2) 
(3.3) 



(3.4) 



where we define orthogonal and longitudinal projections of an arbitrary 3-vector f'' 
with respect to the vector n'' as 



/' = /i + /|l, /|f = ri%/, fi = f-n''njK 
The canonical Hamiltonian of the system is determined by 



(3.5) 



nr 



N 

E 

a=l 



^f^ + e,Ao(t,x,) 
2k^n 



+ ^ / { K (4^^^ + ^o,i)]' + + Ao,iAO'^| d^'x. (3.6) 



Then we get Dirac Hamiltonian which takes into account primary constraints 
and (131): 

X - n^h^'')] d^x, (3.7) 



Hb = h, + J [x{n + n%) + XkiiAnn" + Ao' 



where A, Xk are the Dirac muhiphers. 

The preservation of the constraint ( p.2| ) in time produces the secondary con- 
straint: 

o = {n + n%, H,} = -p- dM + ^ [(Ao'0± - n,y'] 

II 

~ -n\, - p, (3.8) 

where ~ means "weak equahty" in the sense of Dirac and 

N 

p(t,x) = ^e„(5=^(x-x,(t)) (3.9) 

a=l 

is a charge density. 

Next, we consider commutation relations between constraints (P^), (|3.4]), and 



{77 + Wn,, (AttU^ + Ao^^)± - n^h^^} = 0, (3.10) 
{n\i + p, {Ann'' + Ao''=)x - rijh^^} = 0, (3.11) 
{(47ri7'=(t,x) + Ao''=(t,x))x - n^h^'it,^), 
{47rn\t, y) + Ao'\t, y))^ - n^h^\t, y)} 

= _8vr {6'' - n'^n') (n^^^ " y) = ^'^(x - y). (3.12) 
We can check directly, that 

{n\, + p,H,} = 0. (3.13) 

Therefore, the two constraints 

^ = n + Wrii = 0, r = -n\, - p ^ o, (3.14) 

belong to the first class. 

Taking into account (|3.12|) , we come to the conclusion that the constraints 



are of the second class. However, there are two independent second class constraints 
only. Then we can reduce them by means of Dirac bracket: 

{F, G}d = {F, G}- j d^xd'y{F, (47ri7"(t, x) + Ao'"(t, x))x - n,/i*"(t, x)} 
xC„;3(x - y){(47ri7^(t, y) + Ao'^(t, y))x - n,h^''{t, y), G}, a, (3 = 1, 2, (3.15) 
where Cq,/3(x — y) is an inverse matrix to f2°^(x — y): 

j C«^(x - z)W{z - y)dh = 6^J\^ - y). (3.16) 



In the next chapter we shall consider the elimination of the first class constraints 
and the formulation of the Hamiltonian description in the terms of independent 
physical variables. 

Now let us canonically transform the field variables: 



{Ao, n, Ai, Hi) ^ {Ao, S, Ai, £i), 

£ = n + Wrii, Si = Hi, Ai = Ai + AoHi. 



(3.17) 



After the transformation the set of constraints of our system becomes 



S = 0,T = -S\i- p^O, 



0. 



(3.18) 
(3.19) 



Using ( |3.19| ), we can rewrite the canonical Hamiltonian as 



N 



a=l 



ml + k; 
2k„n 



-+eaAo{t,^a) 



in 



(3.20) 



where now kai = Pai - eaAi{t, x^). 



4 Elimination of the gauge degrees of freedom 

Let us consider the first class constraints. We see immediately that Aq, £ are 
a pair of conjugated gauge canonical variables. A second such pair is formed by 
Q = A-^A\i and T, where A = didi, AA'^ = 1, A-^6^{x.) = A-i(x) = -l/(47r|x|). 
Therefore, we can separate the gauge degrees of freedom and gauge- invariant 



ones by means of the Hodge decomposition (see, e.g., [p3], IT 



A=X+diQ, £,=i+^^A-\T + p), (4.1) 

where 

X= (Si - d,A-'d^)A„ i= {51 - dA-'d^)£r (4.2) 

Then we have 

{g(t,x),r(t,y)} = 53(x-y), {A (t,x),i, (t,y)} = A-^)5\^-y). (4.3) 

Since £\i= 0, we can define £^ as follows 

£•= (^L - i^l) ^, e- = vO^f „ = 1, 2. (4.4) 
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Now we have to do a canonical transformation to the new variables ((x^,7raj), 
(0^,60), {Q,r), {Aq,S)), which is generated by the functional 



N 



a=l 



-I 



■ 8 



+ d'A-'{r + p) 



We obtain 



_SF_ 
dF 



61 -6&^ A 



dxi 



d^x. (4.5) 

(4.6) 
(4.7) 



The transverse part of Aii is connected with the new canonical variables as 



Ai= \/4^ (5f - a^A-^a") Oa, a = 1, 2. 



(4.8) 



Now it is convenient to take n = (0, 0, 1) and perform the following canonical 
transformation of the field and particle momenta: 



(4.9) 



The set of constraints for our system in the terms of the new variables can be 
written as follows: 

{E, r, e„ - a„,3) ~ 0. (4.10) 
We reduce the second class constraints by means of the Dirac bracket: 



{F, G}^ = {F, G}- j d'xd^F, e^(t, x) - a^,s(t, x)} 
xC"^/^(x-y){e^(t,y)-a^,3(t,y),G}. 
Here ||C"^(x — y)|| is an inverse matrix to 

W = \\{eait, x) - UaAt, x), ep{t, y) - ap,3{t, y)}|| 



(4.11) 



-25.«^5'(x-y) 



It is given by 



j C-^(x - z)W,p{z - y)dh = S'^S^^ - y). 



C"^(x - y) = -S'^'^Six^ - y^)6{x^ - y^)sgn{x^ - y^), 



(4.12) 
(4.13) 

(4.14) 



because ||C"^(x — y)|| must be antisymmetric. 
Let us introduce the following denotation: 

l5{xy{x')sga{x') ^ ^^^(x). 



(4.15) 
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The final form of the Dirac bracket is 



{F,G}b = {F,G}(^^^^^^) + {F, G}(Ao,£) + {-F, G'}(QT) + -^{F,G}(^a,e) 



6F ^ 6G 



6F 1 6F 



5ea(x) 5ea(x) 6aa{^) ds 6aa{^) 



(4.16) 



where {F, G}(^x,n) denote the standard Poisson bracket in the terms of x and vr. 
Ehmination of constraints (|4.10|) into He leads to the physical Hamiltonian: 



ph 



1 ^ 

-y 

a=l 



1 ^ 



a,fe=l 



Xfol 



-1 1 (a„ - v^|a-v) a (a„ - v^|a-v) rf^^a;, (4.17) 

which generates evolution of an arbitrary function / depending on the gauge-inva- 
riant variables x^, Hai and the terms of the Dirac bracket ( [4.16|) : 



d/ ^ 9/ 
dt dt 



(4.18) 



The gauge-invariant volume element of the constrained field phase space (aa, e") 
is written as 

dr^ph(t) = 7VDetW W Jj5[e,(t,x) -a,,3(t,x)]da,(t,x)de„(t,x), (4.19) 

a=l,2 X 

where 7 is defined as a normalisation constant of the Gauss integral ||^: 

7 j j aa(i,x)L"^(t,x,y)a^(t,y)dW?/^ JJ JJdaa(t,x) 

= Det-i/lL°^(t,x,y)||. (4.20) 

Taking into account (|4.19|) , we can write the volume element of the physical 
phase space of the described system: 



drph(t) = drp,(t)drf,(t). 



(4.21) 



Here 



N 3 



dr^hW=nnd<wd^-w- (4.22) 

a=l i=l 

is the volume element of the particle phase space. 

Now we need check the Liouville theorem: rph(t) = rph(to)- 
It is well known from the classical mechanics, that system evolution in the phase 
space can be described by means of a canonical transformation in the terms of the 
Poisson bracket, which immediately leads to the conservation of the phase space 



11 



volume. This proves the conservation of T^^. In our case the field evolution is 
generated by the Dirac bracket. Nevertheless, the conservation of the volume of the 
constrained field phase space can be proved and we shall demonstrate its possible 
proof for the considered Dirac bracket ( [4.16| ) in Appendix B. 
The volume element of the full phase space is 

drMi(t) = drph(t)drg(t), (4.23) 

where drg(t) is the gauge volume element of the phase space: 

Yl S[S{t, x)](5[r(t, x)]dAo(t, x)d^(t, x)dQ(t, x)dr(t, x). (4.24) 

X 

However, we turn off drg(t) from the further description, because all thermodynam- 
ical characteristics of the system do not depend on dynamics of the gauge degrees 
of freedom (see [0). 



5 Statistical mechanics 

Now let us imagine that the physical initial data xl^ito), vrai(to)) CLa{to,y^), eQ,(to,x) 
are not precisely known. Hence, we can introduce a probability density ^(to) = 
g{to, x^(to), 7raj(to), CLa(to,yi), (to ) x)) for having various initial states. This func- 
tion satisfies the condition 

J f?(to)drph(to) = 1. (5.1) 
Then average value of a general dynamical variable / is defined by 

7W = j fit, <it), ^ai{t), a,(t,x), e„(t,x))^(t)drph(t). (5.2) 

Since to has been randomly selected, we come to the relation 

1 = J Q{to)dTpM = J f?(^)drph(t). (5.3) 

We have already seen that drph(to) = drph(t), so one immediately obtains the 
Liouville equation 

= '-§='-j^+m,HMt)W (5.4) 

Taking into account ( p.4| ), we interpret g as an integral of motion. 

Let us consider the canonical Gibbs ensemble in equilibrium. In this case we 
have 

^(t) = Ce-^^i"^(*), (5.5) 

where f3 = 1/kT and C is a normalisation constant. 
Partition function can be found as 

-/3Hph(i) 

drph(t). (5.6) 
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We shall find below the value of Z. 

Let us first rewrite the physical Hamiltonian as follows 



ph 



1 ^ 

4e 



a=l 



1 ^ 



|Xa - Xb 



+ 



27r 1 ^p(t,x)A-i(x-y)|p(t,y)rf='xA 

27r 1 (|p(t,x)+j°(t,x)) G(t,x-y) (|p(t, y) + J°(t, y)) rfWy 



where 



a„(t, x) = a„(t, x) - x - y) (^^p(t, y) + y)^ dS, (5.8) 

(5.9) 



J°(t,x) = -X:e.^5^(x-x.(t)) 
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G(t,x) = (5^(x). 



(5.10) 



Since ifph does not depend on field momenta e°, after an integration over in 
( p.6| ) we obtain the expression for the partition function: 



1 ^ 



|Xa - Xfol 



1 /■ 1 / 

Z = 7v/DitW— / exp(-/5) -- 5^ vr,3 + 

■ L a=l ^ 

-27r 1 ^p(t,x)A"i(x-y)^p(t,y)ciW2/ 

+2vr 1 (^|p(t,x)+j°(t,x)^ G(t,x-y) (^|p(t, y) + jO(t, y)^ ^3^^^^ 
-i y S„(t,x) ^A + 47r^-^53(x-x,(t))'j a,(t,x)rf3 



nn^ n n^^-c-) 



a=l i=l 



o=l,2 X 



Here we can replace daf^ by da^. Taking into account 
yields 



(5.11) 

the integration over da 



Z = Z^ 



N 3 



M/nn 



a=l 1=1 



2n 



exp(-/3) 



a=l ^ 



TTaS 
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^ T^, — Xb| J C3 O3 



a, 6=1 
+27T 



5, 



p(t, x) + f^{t, x) ) G{t, X - y) ( ^p(t, y) + £{t, y) ) rf^x^ 
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X 



ei 



6^^ - y) + 47rA-i(x - y) V ^6%y - x„) 



where 



VDetW 



VDetll -/55"/^A53(x-y)| 
represents the free field partition function. 



(5.12) 



(5.13) 



6 Conclusions 

In this paper we have once more demonstrated the usefulness of the various forms of 
relativistic dynamics in treating the relativistic particle systems. Constraint analy- 
sis of the considered problem shows significantly different structures of Hamiltonian 
description of charged particles, interacting by means of an electromagnetic field, in 
the space-like and isotropic forms of dynamics. Specifically, the use of front form of 
relativistic dynamics allows us to exclude the electromagnetic field variables from 
classical partition function of the system of charged particles. It is obvious that 
such an exclusion must depend on the particular boundary conditions for the field 
variables, but at present the proper sense of these conditions remains unclear. The 
obtained representation for the partition function contains highly nonlocal expres- 
sions, and their further analysis consist of a complicated task. The consideration of 
the various approximation schemes at this step seems to be inevitable. 

It should be noted that exclusion of the field variables transform the problem into 



the domain of relativistic direct interaction theory (see, e.g., p^, ^). An application 
of such a theory to the consistent formulation of the relativistic statistical mechanics 
is just at the beginning. 

On the other hand, the established Liouville equation may be used in nonequi- 
librium situation as well. It serves an useful starting point for deriving various new 
forms of kinetic equations for charged particle system. 



Acknowledgment 

We are greatly indebted to Yu. Yaremko and V. Shpytko for many stimulating 
discussions. Numerous helpful conversation with A. Duviryak and his support in 
the performance of this work are especially acknowledged. 



14 



Appendix A 

Relativistic free particle in an arbitrary form of 
dynamics 



The Lagrangian of relativistic free particle in a given form of dynamics (|2.5|) is 

(A.l) 



-m\/ {DifY — = —niT ^. 



The canonical momentum and Hamiltonian are given by 

Pi = mT{gijV^ - ipiipt), (A. 2) 

H = mViftDip, (A.3) 

where the matrix gij = Sij — (fi^pj has been introduced. The determinant of the 
matrix is 

g=\\g^^\\ = l-<p^, (A.4) 

and (? > as result of the condition ( |2.10| ). 



Consider firstly the case g > 0. Then the matrix gij has an inverse, gij 
6ij + g~^(pi(pj and from (|A.2|) it follows 



QijPiPj = {^^9 Vt - 1) 



Using (|2.8| ), we find 



Vi = B^^ipt [pi + g~^ipi{B + Pjip^)] , 
Combining these results with (|A.3|) , one gets 



H = iptg 



-1 



g{m'^ + p2) + {piif'y + pi(p' 



(A.5) 

(A.6) 

(A.7) 
(A.8) 

(A.9) 



The case of isotropic forms of mechanics can be treated by taking the limit g +0 
and using that, in view of ( |A.2| ), 

(A.IO) 



Piip' = mT [gviip' + ifitig - 1)] 
so that piip^ <0 as g ^ +0. It gives 

1 + p2 

H = --ft — • 

2 pi(f^ 



(A.ll) 
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Appendix B 

Proof of the Liouville theorem for a given Dirac 
bracket 

Let us show that evolution of the field variables as the generalized canonical transfor- 
mation in the terms of the Dirac bracket conserves the volume of the constrained 
field phase space. 

We first consider transformation generated by some functional G with an arbi- 
trary parameter S,'- 



Oia — da + ^{tta, G}d — O-c 



i5G_ _ i 1 5G 



^,1 ri f6G 



(B.l^ 



We immediately see that — a^.s = Cq, — aa,3- So, the constancy condition of 
Fpjj leads to the relation: 

/ n n '^[^"(^' ^) ^ "a,3(i, x)]daa(t, x)dea(t, x) 

a=l,2 X 

= [ Ji^,0) n n'^t^"^^'^) -^"'3(t,x)]da„(t,x)de,(t,x). (B.2) 



a=l,2 X 

Here Jacobian J(^, t]) is defined by 



J(e,r/) =Det 



(5(a„(t,x;0ea(^,x;0 



If ?7 = 0, then we have 



J(^, 0) = Det 



(5(aa(t,x;0ea(t,x;0 



5(a^(t,y)e^(t,y)) 

It is evident that J(0, 0) = 1. Let us compute (dJ(^, 0)/d,^)L q. We get 



(B.3) 



(B.4) 



dJ(e,o) 


= Tr 


59gaQ(t,x;0 


+ Tr 








5=0 








5=0 



= Tr 
-Tr 



6^G 



2 6ap{t,y)6ea{t,yi) 
1 



2 5e^(t,y)5a„(t,x) 



Tr 
f Tr 



1 1 



5^G 



2d3 6ai3{t,y)6aa{t,x.) 
5^G 



-(93 — 

2 (5e/3(t,y)(5e„(t,x) 



(B.5) 



It is obvious that the first and the third terms vanish. One of possible ways demon- 
strating a cancellation of the second and fourth terms is based on commutation of 
the Dirac bracket action and replacement of 3 by (or Cq, by a^^^ in view of 
- ao,3 = 0). 
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Therefore, we have 

dJ(e,o) 



Since J(e,0) = 7(^,6)^(6,0), then 



= 0. (B.6) 



dJ(e,o) 



^ dJ(66) 



J(6,0) = 0. (B.7) 



Thus, J(6 0) = 1 for all 6 

Now if we take ^ = t and G = H^^, we come to conclusion 

rph(t) = rph(to), (B.8) 

namely, the time evolution preserves the phase space volume Fph- 
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